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Abstract 

We study the persistent current of noninteracting electrons subject to a pointlike magnetic flux in 
the simply connected chaotic Robnik-Berry quantum billiard, and also in an annular analog thereof. 
For the simply connected billiard we find a large diamagnetic contribution to the persistent current 
at small flux, which is independent of the flux and is proportional to the number of electrons (or 
equivalently the density since we keep the area fixed). The size of this diamagnetic contribution is 
much larger than mesoscopic fluctuations in the persistent current in the simply connected billiard, 
and can ultimately be traced to the response of the angular momentum I = levels (neglected in 
semiclassical expansions) on the unit disk to a pointlike flux at its center. The same behavior is 
observed for the annular billiard when the inner radius is much smaller than the outer one, while 
the usual fluctuating persistent current and Anderson-like localization due to boundary scattering 
are seen when the annulus tends to a one-dimensional ring. We explore the conditions for the 
observability of this phenomenon. 

PACS numbers: 73.23.Ra, 73.23.-b, 73.43. Qt, 75.75.+a 
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I. INTRODUCTION 



A resistanceless flow of electrons can occur in mesoscopic systems if the linear size L is 
less than the phase coherence length L$. The simplest example of this is a one-dimensional 
metallic ring threaded by a magnetic flux $. The thermodynamic relation 

dF 

defines the persistent current in MKS units. At zero temperature, which we will focus on, 
the free energy F can be replaced by the total ground state energy E. 

Persistent currents were first predicted to occur in superconducting rings 1 ' 2,3 . It was later 
realized that persistent currents exist in normal metallic rings as well 4,5 . The phenomenon is 
understood most easily at zero temperature for a ring of noninteracting electrons, where the 
electronic wavefunction extends coherently over the whole ring. If the ring is threaded by 
a solenoidal flux, all physical properties are periodic in applied magnetic flux with a period 
of the flux quantum $ — h/e. A nonzero flux splits the degeneracy between clockwise 
and anticlockwise moving electrons. Upon filling the energy states with electrons, one finds 
ground states which have net orbital angular momentum, and net persistent current. Much 
experimental work has been carried out on ensembles of rings/quantum dots 6,7 in a flux 
as well as on single metallio^ 1 ^^ 2 - or semiconductor quantum dots/rings^ 3 - 1 ^ 4 - 1 ^ 1 ^. The 
subject has a long theoretical history as well^^^^^^^^ 4 ^ 1 ^. 

In this paper we investigate the persistent current of noninteracting electrons in quan- 
tum billiards subject to a point flux. Related semiclassical calculations have been car- 
ried out in the past for regular (integrable in the absence of fiux^^^ and chaotic 
billiards^ 7 -* 2 ^^. Numerics have previously been performed on these systems as well— >2&. 
We carry out calculations on the simply connected chaotic Robnik-Berry billiard 3 ^^ 7 - 1 ^'^' 4 ^ 
obtained by deforming the boundary of the integrable disk, and on an annular analog which 
we call the Robnik-Berry annulus. The ratio of the inner r to the outer radius R of the 
annulus (£ = r/R) plays an important role in our analysis, and allows us to go continuously 
between the simply connected chaotic two-dimensional billiard and a (effectively disordered) 
quasi-one-dimensional ring. 

Our main result is that there is a large diamagnetic and flux-independent contribution to 
the persistent current for |$| $o in the simply connected billiard which is proportional to 
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the number of particles, and overwhelms the mesoscopic fluctuations which have been the 
focus of previous work— t M i ^ i Si i ^ ^ I^^^o^i^ ^i M i ^. This arises from angular momentum 
I = states in the integrable disk, which respond diamagnetically, with energy increasing 
linearly with the point flux, for small flux. This behavior is robust under the deformation of 
the boundary which makes the dynamics chaotic. As £ increases from zero, this contribution 
to the persistent current persists for typical $/$o — 1> but smoothly decreases in magnitude 
and becomes negligible for £ — > 1. The precise £ at which the diamagnetic contribution to 
the persistent current becomes equal to the typical fluctuating paramagnetic contribution 
depends on the electron density. For £ 7^ and very tiny flux the diamagnetic contribution 
to the persistent current varies linearly with <&/<&o ( see below) . 

This diamagnetic contribution seems to have been missed in previous work, to the best of 
our knowledge. The reason is that the semiclassical approximation becomes asymptotically 
exact as the energy tends to infinity, and in this limit, the spectral density of I = states 
vanishes. Thus, I = states are explicitly disregarded^^ 1 ^ 1 ^ i n the semiclassical approach, 
since they do not enclose flux. It has been noted in the past that diffraction effects necessitate 
an inclusion of I = states in the sum over periodic orbits on the integrable disk—, but the 
connection to persistent currents was not made. 

It should be emphasized that since the total persistent current is a sum over the con- 
tributions of all levels, the diamagnetic contribution we uncover exists even at very large 
energies, where the levels at the Fermi energy are well approximated by semiclassics. 

The robustness of the diamagnetic contribution to the persistent current under deforma- 
tion can be understood as follows: In the chaotic billiard, each state at a particular energy 
is roughly a linear combination of states of the disk within a Thouless energy (E? — hvp/L, 
where L is the linear size of the billiard) of its energy. When the Fermi energy Ep greatly ex- 
ceeds E T , the contribution of the occupied states does not change much when the boundary 
is deformed and chaos is introduced. 

This behavior appears similar to, but is different from Landau diamagnetism^ in a fi- 
nite system, which is a response to a uniform magnetic field. The primary difference is 
that the orbital magnetization (proportional to the persistent current) in Landau diamag- 
netism is proportional to the field itself (because the energy goes quadratically with the 
field strength), whereas the effect we describe is independent of the flux for small flux in the 
simply connected Robnik-Berry billiard (because the energy goes linearly with the flux). In 
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the Robnik-Berry annulus with £ — > 0, the energy rises quadratically with the flux for very 
tiny flux $ <C $o/logA£ -2 , but crosses over to the linear behavior characteristic of the 
simply connected system for larger <3>. Since the flux is pointlike, and in the annular case 
nonzero only where the electron wavefunctions vanish, the entire effect is due to Aharanov- 
Bohm quantum interference. Since the effect is primarily caused by levels deep below Ep, 
experimental detection is feasible only through the total magnetization, and not by conduc- 
tance fluctuations which are sensitive to the levels within the Thouless shell (lying within 
Et of Ep). Previous samples have been subjected to a uniform field rather than a point 
flux-^^iit^^i, anc i anyway the ring samples have £ too large for this effect to be seen. 
However, we believe that experiments can be designed to observe this effect. 

The plan of the paper is as follows: In Section II we describe the method we use to cal- 
culate the spectrum, and present analytical expressions and numerical results for persistent 
current in the disk and simply connected Robnik-Berry billiards. In Section III we gener- 
alize the method to the annulus and present our results. Conclusions and implications are 
presented in Section IV. 



II. THE SIMPLY CONNECTED ROBNIK-BERRY BILLIARD 

We begin by briefly describing the procedure to obtain the energy levels within the 
billiard, which leads to the persistent current: 

k 

We work with the Robnik-Berry billiard^ 1 ^, which is obtained from the unit disk by 
conformal transformation. The original problem of finding energy levels of electron in the 
domain with complicated boundaries is reduced to a problem where the electron moves in 
the unit disk in a fictitious potential introduced by the following conformal transformation; 

. , z + bz 2 + ce l5 z 3 
w(z) = . (3) 
V ; VI + 2b 2 + 3c 2 V ; 

where w = u+iv represents the coordinates in the laboratory coordinate system, and z = x+ 
iy are the conformally transformed coordinates (details are in Appendix A). The parameters 
b, c, and 5 control the shape of the original billiard, and for the values we use, the classical 
dynamics is mixed, but largely chaotic. It is also straightforward to introduce a point flux 
which penetrates the center of the unit disk— ^§j22ii£ (after the conformal transformation). 
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FIG. 1: Ground state energy Eg in units of 10 4 x 2r ^s i an d persistent current I in units of 47r ^ R 2 
as a function of dimensionless flux for the regular disk (panels a,b) and the simply connected 
chaotic billiard (panels c,d). The results are for 200 particles. 

We find 600 energy levels for regular and chaotic billiards for different values of parameter 
a that controls magnetic flux coming through billiards. Only the lowest 200 levels are actu- 
ally used in further calculations, since the higher levels become increasingly inaccurate^ 4 ^. 
The persistent current is obtained as a numerical derivative of the ground state energy for 
a given number of electrons. 

For the unit disk billiard the ground state energy Eq has a non zero slope as a — > 0. 
Thus there is a persistent current in the system for arbitrarily small magnetic flux (see Fig. 
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Qualitatively this behavior can be understood as follows. In the absence of magnetic field 
energy levels corresponding to orbital quantum numbers ±Z, are degenerate. A nonzero $ 
lifts the degeneracy and for small a the two ±/ levels have slopes that are equal in magnitude 
and opposite in sign. Thus, as long as both are occupied, these levels do not contribute to 
the net persistent current /. The only nonzero contribution comes from levels with I = 0. 

For the unit disk the expression for the persistent current can be derived analytically for 
small values of magnetic flux. At zero temperature the persistent current due to kth level 
is Ik = —d€k/d& (e*; is a dimensionless energy, and Ik is persistent current divided by the 
energy unit h 2 /2mR 2 ; see appendix |A] for notations). For the unit disk, the energy levels 
are found from the quantization condition: 

J\l-a\{l\l-a\,n) = 0, 6 k = - a\). (4) 

Then from Eq. ([21) , persistent current caused by kth level is: 

h = -jln(\h-a\) . (5) 

To find dj/da we differentiate Eq. (TJJ: 

dJjg) = dJvin) dv , dJ v (7) 9j = Q , Q s 
da dv da d'y da 

For I = levels, v — \l — a\ — a. In a — > limit, for the derivatives of Bessel function 



one gets: 



dMi) 



dv 
dUl) 



c?7 



= > ( 7 ), 

,=o (7) 

= -Hi). 

Combining Eqs. (jSJ) and ((Zj) and using relation that for 7 ^> 1 Bessel function iVo( 7 ) ~ 
Ji( 7 ) (this approximation works well already for the first root of Eq. (jlj)), we find §^| a =o = f , 
which leads to: 

/ = -xE^(°)' ( 8 ) 

n 

where summation is over the levels with orbital quantum number I = 0. 

For large argument values (which is the same as large energies) the quantization condition 
(BJ for the unit disk becomes cos(7„ — na/2 — 7r/4) = 0, with roots: 

7„ = 7ra/2 + 7r/4 + 7r(2n + l)/2. (9) 



With the energy being measured in h 2 /2mR 2 units, the Fermi wave vector is kp = ^max ~ 
Kn max , where n max denotes the largest / = level. With disk area equal to 7r (R — 1), the 
number of particles in the system is iV = {im max /2) 2 . This allows us to find the dependence 
of the persistent current on number of particles in the system in a —>■ limit. 



where we neglected a subleading term proportional to n max . We remind the reader that the 
physical persistent current is the expression in formula fflUj) divided by energy unit h 2 /2mR 2 . 

In Fig. [2jthe persistent current I is plotted against the number of particles iV for magnetic 
flux a = 0.01. The behavior of the current is consistent with Eq. ([TO]) . That is, for small 
magnetic flux it is proportional to 2N. For the regular disk (Fig. [2k) the persistent current 
is a set of consecutive steps. Each step appears when the next I = level is added to 
the system. The length of the steps is equal to the number of I ^ levels between two 
adjacent levels with zero orbital quantum number. As one particle is added to the I ^ 
level, it results in persistent current jump. The next level has opposite slope, and once it is 
occupied, cancels the contribution of the previous I ^ level to the net persistent current. 
This explains the noise above each step in Fig. Eh- In addition, each step has a small 
inclination which is due to the fact that the I ^ levels do not cancel each other exactly 
when $ 7^ 0. For larger magnetic flux the steps become more inclined. 

To see that levels with I ^ do not contribute to persistent current at weak magnetic 
flux we simply note that derivative of 7: 



is an odd function of /. In the a — > limit, the root 7„(^) and the derivatives of J v (7 n ) in 
Eq. ffTTj) are even functions of I, and dv/da is odd. As a result, the whole expression is odd 
function of I, which proves the cancellation of ±Z levels in Eq. ([5]). 

For the chaotic simply connected Robnik-Berry billiard each eigenstate is a superposition 
of all /-states of regular disk (see Eq. flA5j) ). mostly within a Thouless shell of its energy. 
Assuming that states with ±Z enter this superposition with equal probability over the en- 
semble due to the chaotic nature of motion, one can conclude that the ensemble-averaged 
contribution of these levels to the net current is zero. However, as seen in Fig. [2b, the meso- 
scopic fluctuations due to the I 7^ levels are overwhelmed by the diamagnetic contribution 
linear in iV for small a. 
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FIG. 2: Persistent current / vs. number of particles N for regular disk (a) and chaotic disk (b) for 
the value of reduced magnetic flux a = 0.01. 

III. THE ROBNIK-BERRY ANNULUS 



Now we turn our attention to annular billiard. Here there is an additional parameter £, 
which is the ratio of the inner radius r to the outer radius R of the regular annulus in the 
(conformally transformed) z plane. By varying £ we are able to smoothly go from the simply 
connected Robnik-Berry billiard to an (effectively) disordered ring in the limit £ — > 1~. 

First consider the disk limit £ — > 0. We can derive an analytical expression for I when 
£ is small enough that £7„ <C 1, and 7„ > 1. For a regular annulus with _R = 1 energy 
quantization follows from the Dirichlet boundary condition: 



Uln)N u ( ln i) - J„( 7 nOW„(7n) = 0. 



(12) 



Here n numerates root at fixed angular momentum v. We use the large and small argument 
expansions for Bessel functions to obtain for the / = levels: 



1 



not it 
cot( 7n - — - -) 



rfl + a) V 2 



2 2 cot(aTr) { ln€\ a r(a) f ^ 

2 



(13) 



r(l + a) V 2 j it 

We express the roots for the annulus as a small deviation from the roots for the disk, 
which we denote 7^; 7„ = 7^ + #7„ with 7^ = air /2 + tt / '4 + n(2n + l)/2. Approximating 
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cot(a7r) by l/(a7r), for small values of 8^f n we find: 

an 



1 + coth a In 



7n£ 



(14) 



In Eq. dHJ), for small a, 7 n under logarithm can be safely replaced by its value for the 
disk jn ■ Then roots for the annulus are: 



7T 7T . 7T , 

7„. = - + 2"( 2ri + 1) — 2 a c ' a 



(15) 



One can now take various limits of Eq. (|15[) . To recover Eq. for the disk roots, we keep 
magnetic flux a fixed and take the limit £ — > 0. As one can see from Eq. (fT5]) , convergence 
to the disk limit is slow due to the logarithm, and occurs only for a>l/ log (7n £/2). 

Another limit of interest is to keep £ fixed and obtain behavior of roots 7 n for small a. For 
small a C 1/ log (7^ £/2), the roots 7„ with / = vanish quadratically with a. Expanding 
the coth function in Eq. ( ITBT) . we obtain: 



7T 7T , 

7 «= 4 + 2 (2n + 1) " 
which leads to the persistent current: 

2irea 



71 



a 



- 1 + — In' 



7^ 



In 



3A 



E 



- + -(2n + l' 
4 2 V 



In 



7^ 



-1 7n s 



Rough estimation of this sum with help of the Euler-MacLaurin formula gives: 



7i 2 ea 9 

I f~ 



1 Tlmax^'^ 
Ti" \W 

3h max 2^i ' 



(16) 



(17) 



(18) 



where we kept only terms proportional to n max and e = 2.71828... inside the logarithm 
denotes Euler's number and not the electronic charge. 

Using the relation N = (iTn max /2) 2 (for small values of £ the density of states for the 
annulus and the disk are practically the same), the persistent current becomes for small 
a«l/log( 7 ^/2): 



a 



In 



/(<*) = -I2N. 

h 



(19) 



To approach the limit of a one-dimensional ring, where 7 n ^> 1 and 7„£ ^> 1, we return 
to quantization condition (fl2l) and use the following large argument expansion for Bessel 
functions: 



We use formulas (!20|) and quantization condition (jT2l) to get a new equation for roots: 

v 2 - 1 /4 

sin(7 n a) - cos(7„cx) — — — a = 0, (21) 

where we ignore the term proportional 1/7^. The quantity a = 1 — £, is assumed to be 
much less than unity. For sufficiently small a one can drop the second term in Eq. (l2Tj) and 
get 7 n <j = 7m. To find corrections to this expression we assume that 7 n <7 = irn + rj with 
rj = - 2 ~^ 4 cr 2 <C 1 and plug it in Eq. (|2"T|) to obtain the solutions of quantization condition 
(J2Q are: 

7m z/ 2 - 1/4 ., . . 

7» = — + — <t, i/=i-a. 22 
a 2irn 

The energy spectrum for the annulus in this limit is: 

^ = 7^(^) 2 + (^ 2 -l/4). (23) 

The first term in Eq. (1231) denotes the radial kinetic energy and diverges in a — > limit. 
This divergence can be absorbed into the chemical potential for the n = 1 radial state. The 
difference between energy levels with radial quantum numbers n is of the order n(ir/a) 2 . 
For cr —> => £ — > 1 one can assume that all the levels of interest have the radial quantum 
number n — 1, and are labelled only by orbital quantum number /. Since our diamagnetic 
persistent current arises from a large number oc of / = levels, it is clear that it vanishes 
in the limit of a ring. 

It is straightforward to show that for a regular annulus the contributions of ±/ levels also 
cancel each other for small values of a. However, levels with I = have zero slope when 
a — > 0. To show this one takes the derivative of quantization condition f|T2|) : 



+ 



da 



UMMln) + Ulnt)N v { ln ) - J v {ln)N v { ln i)-iJ v { ln )N v { ln £) 



0, (24) 



where A v (z) = dA v (z)/du, and A' u (z) = dA„(z)/dz. When a — > 0, derivatives of 
Bessel functions become J^(z) = ttN (z)/2, N v {z) = —ttJ (z)/2, J' v {z) = —Ji(z), and 
N' v {z) = —N\(z). Then all terms outside square brackets in Eq. ( 1241) cancel each other. The 
expression inside brackets in general has a non-zero value, which means djn/da = 0. 

In Fig. [3] the persistent current in the annular billiard is depicted for different values of 
the aspect ratio £. To facilitate the comparison between different values of £, we keep the 
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FIG. 3: Persistent current I vs. reduced magnetic flux a for several values of £ for N = 200 
particles. Figure (a) represents current for regular annulus normalized to the same density of 
states (same area) for different values of £. The current for the chaotic annulus is depicted in 
Figure (b). 

area of the annulus the same, thus keeping the average density of states the same. For a 
regular annulus (Fig. (3^) for small values of flux, the current is a linear function of a. As £ 
gets smaller, the diamagnetic contribution to the persistent current increases. This behavior 
is consistent with Eq. (fT9l) that shows linear dependence on a and slow growth as £ — > 0. 

In the regular annulus, for £ close to unity the behavior of the persistent current is close, 
but not identical, to that of a ID ring. Even for £ = 0.9 there exist several states with 
1 = 0, which means that our billiard is not purely a ID ring. The effect of these states on 
the persistent current is not entirely trivial. For a fixed number of particles in the system, 
as £ changes, the number of I = levels also changes. As the next / = level is added (or 
expelled), the current experiences a jump. The magnitude of this jump is large enough for 
small a to make current to be positive. For larger a the current remains diamagnetic. 

In the distorted annulus (Fig. [3)d) the persistent current is a linear function of a for small 
a. For larger magnetic flux one observes nonlinear behavior that can be attributed to level 
repulsion in the chaotic billiard. 

The dependence of the persistent current in the annulus on the number of particles iV at 
fixed a is similar to that in the simply connected billiard. At small £ the persistent current 
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FIG. 4: Persistent current I in distorted annulus vs. number of particles N for several values of £. 
magnetic flux a = 0.17. 
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FIG. 5: Persistent current I vs. number of particles N for distorted annulus. Parameters b, c, and 
5 control the shape of billiard. 

in the regular annulus is a staircase-like function. For the distorted annulus the numerics 
are scattered around a straight line (see Fig. Hj). 

For larger £ the magnitude of diamagnetic contribution to the persistent current decreases, 
and the numerics are dominated by mesoscopic fluctuations. When £ — > 1, the persistent 
current becomes negligible (see Fig. [5^) for low occupations. We believe this is a manifes- 
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tation of Anderson localization due to the boundary scattering. At high energies, when the 
localization length exceeds the circumference of the annulus, extended states reappear and 
can carry persistent current. In Fig. Owe plot the persistent current for 2 different sets of 
parameters controlling the shape of annulus. For a large distortion (Fig. the current is 
nonzero only for high energy states beyond N = 110. In Fig. \5jp the parameters b, c, and S 
are chosen to make the annulus less distorted, and we see that the threshold for extended 
states moves to lower energy (about iV = 40). 

IV. CONCLUSIONS, CAVEATS, AND OPEN QUESTIONS 

We have investigated the behavior of chaotic simply connected and annular billiards 
penetrated by a pointlike flux. The annular billiards are characterized by a dimensionless 
aspect ratio £ = r/R, the ratio of the inner (r) to the outer radius (R). Note that in the 
annular billiards, the flux exists in a region where the electrons cannot penetrate, and the 
effects of the flux on the electrons are purely Aharanov-Bohm quantum intereference effects. 
Our main result is that there is a systematic diamagnetic contribution to the persistent 
current which can be traced back to the flux response of the I = levels of a regular unit disk 
(or annulus). Even though the number of such I = levels is submacroscopic (oc \fN, where 
N is the number of electrons), the contribution to the persistent current due to these levels is 
proportional to N and is independent of the flux for small flux in simply connected billiards, 
and can overwhelm the fluctuating mesoscopic contribution^^^^^^^^^^ 1 ^ 1 ^ 1 ^ from 
the states in the Thouless shell (\E — Ep\ < Ej-). This effect is quite distinct from Landau 
diamagnetism^i. 

The diamagnetic contribution to the persistent current from I = states seems to have 
been missed in previous work using the semiclassical sum over periodic orbits^^ 1 ^ 1 ^. This 
is understandable, since the semiclassical approach becomes exact only as E — > oo, and in 
this limit, the I = states have vanishing spectral density pi = o(E) ~ 1/yE. However, we 
emphasize that the total persistent current contains the sum over all levels, and will indeed 
behave diamagnetically at small flux (in the simply connected billiard) as we have described. 

For very tiny £, the annular Robnik-Berry billiard behaves much like the simply connected 
one for most values of the dimensionless flux a = $/$o ^ 1/ log A£~ 2 , with a diamagnetic 
contribution to the persistent current which is proportional to the electron density. However, 
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convergence to the £ = limit is logarithmically slow, and the limits a — > and £ — > do 
not commute. As the aspect ratio £ increases, and the annulus tends to a one-dimensional 
ring, this effect diminishes to zero. For £ close to 1, we also see Anderson localization 
in the distorted annular billiards, wherein the persistent currents are negligible below a 
certain energy (presumably because the localization length for these levels is smaller than 
the circumference), and become nonzero only beyond a threshold energy. 

While we can obtain analytical estimates for the limits £ — > and £ — > 1, it is difficult 
to make analytical progress for generic values of £ (not close to or 1). However, one can 
easily verify from the asymptotic expansions that for generic £ the diamagnetic contribution 
to the persistent current for $ -C $ g° es as 

(25) 

where r, R, are the inner and outer radii respectively. This should be compared to the 
typical fluctuating persistent current for interacting particles 1 ^^ 1 ^ which behaves as 




E T n 2 N 



l flue 



(26) 



$ mR% y R(R - r) 

It can be seen that the ratio of the systematic diamagnetic persistent contribution to the 
fluctuating contribution is roughly 

l J **l ~ ( R ~ r ) $ (27 ) 
\Ifiuc\ r $ 

Previous ring sample:^ 4 ^ 4 ^^^ 4 ^ have (R— r) <r. They are also subject to a uniform 
magnetic field rather than a point flux. Despite this, a systematic diamagnetic contribution 
at low flux has been detected in recent experiments^ 1 ^. However, the experiments are 
carried out at finite frequency, and the effects of attractive pair interactions^ 1 ^ (see below) 
or nonequilibrium noise 44 cannot be ruled out. 

In order to detect this effect unambiguously, one must work with a material which has 
no superconductivity at any temperature, to rule out attractive pair interactions. It is also 
clear that needs to be made as large as possible in order to render this effect easily 
observable. Care must be taken that there is no magnetic flux in the region where the 
electron wavefunctions are nonzero in order to maintain the pure Aharanov-Bohm quantum 
interference nature of this effect. 
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Let us now mention some caveats about our work. We have taken only a few (« 200) 
levels into account, whereas most experimental samples have a hugely greater number of 
levels. However, the physics of the diamagnetic contribution to the persistent current for a 
particular level concerns only whether that level has I = or not, and is independent of its 
relative position in the spectrum. We expect our conclusions to hold for arbitrary densities. 

We have considered a pointlike flux, which is unachievable in practice. For the annular 
billiard, all one needs to ensure is that the flux is nonzero only in the central hole of the 
annulus, and is zero in regions where the electron density is nonzero. By gauge invariance, 
such a situation will be equivalent to the one we study. 

We have also ignored the effect of interelectron interactions. For weak repulsive 
interactions^ 1 ^ 1 ^ 1 ^, we expect interactions to modify the effect only slightly, because it 
comes primarily from occupied levels deep within the Fermi sea, which are Pauli-blocked 
from responding to the interactions. However, for strong repulsive interactions^^*^, sig- 
nificant corrections to the persistent current^ from electrons in the Thouless shell cannot 
be ruled out. If the interactions are weak but attractiv^^^^, the low-energy fluctuations 
of Cooper pairs become very important^ 1 ^, and can produce additional large diamagnetic 
contributions at low fields. 

Similarly, though we have concentrated on the zero-temperature behavior, we expect this 
effect to persist to quite high temperatures, since most of the I = levels involved lie deep 
within the Fermi sea. 

Finally, it would be interesting to investigate the effects of static disorder within a chaotic 
billiard, which would induce the system to cross over from a ballistic/chaotic to a disordered 
(diffusive) system. We hope to address this and other issues in future work. 
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APPENDIX A: NUMERICS FOR ENERGY LEVELS 



The idea of this method is as follows^ 1 ^ 1 ^ 1 ^. In the original (uv) domain the Schrodinger 
equation is: 

1 9 

— (-ifiV - qA(u,v)) = EV(u,v), q = -e<0. (Al) 

To keep dynamics of electron unchanged, it is assumed that magnetic field exists only at 
the origin of (uv) plane inside the billiard. This requires that vector potential satisfies the 
condition V x A(r) = n<M(r), where n is a unit vector perpendicular to the plane of the 
billiard. 

The billiard is threaded by single magnetic flux tube. The strength of the flux is $ = a$0) 
where $o — h/e is a magnetic flux quantum. 
If the vector potential has the form: 

A ^ = ^"{%-- 9 /u- )- ! =\^< < A2) 

then with the help of conformal transformation: 

, , z + bz 2 + ce i5 z 3 . . 

w(z) = — , w = u + iv , z = x + iv. (A3) 

1 ; VI + 2b 2 + 3c 2 K J 



the Schrodinger equation in polar coordinates of (xy) plane becomes: 

—d e ^(r,9) - - 



V^(r, 9) - —d e ^(r, 9) - — *(r, 9) + e w (re ie ) *(r, 9) = 0. (A4) 



Here the energy e is measured in units of h 2 /2mR 2 , and the distance is in units of R, where 
R is the radius of the disk in (xy) plane. Also, the coefficients b, c, and S in Eq. ( 1A3j) are real 
parameters selected in the way so that > for all values of z inside the disk in (xy) 

plane. The transformation w(z) is a cubic polynomial normalized to preserve the area of 
the billiard and leave the density of states invariant. Equation flA4j) should be accompanied 
by Dirichlet boundary condition. 

To find the energy spectrum, one expands the \?(r, 9) function in Eq. ( 1A4I) in terms of 
the eigenstates <pi t n(r, 9) of free electron(w = 0) inside the round billiard (R = 1): 

oo (p) 

%(r,0) = N p Y t -*-<l> j (r,e). (A5) 
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Compound index j = (u, n) numerates levels in ascending order. Normalized function 

(f>in(r,Q) is: 



J\l-a\(l\l-a\,nr)e l 



(A6) 



V^j'\l- a \{l\l-a\,n) 

Function J v {r) is the Bessel function of the first kind, 7^ is the nth root of Jv(r), and / is 
an orbital quantum number. The coefficients of expansion in Eq. ( 1A5I) are chosen that way 
for further convenience. 

Plugging expansion (1A5|) into Eq. flA4p . after simplification one gets the matrix equation 
for eigenvalue problem: 

(A7) 



M l3 cf 



» 



where matrix M is: 



Ma 



+ S iaj . 2 6ce- lS I^ + <U-i(4My + Ubce^l, 



(i) 



-i5r(3)> 



+ KiA^if + 18c 2 4 4) ) + S liilj+1 (4b4 ] + 12bce i5 I®) 
+ 5 hih+2 Qce l5 lg ) 



(A8) 



/(l + 26 / + 3cf). 



The integrals iff have the form: 



j(h) = So drr^JvAlirfJujiw) 
13 liljJLtidJlAlj) 



(A9) 



Along with the simply connected domain (irregular disk) we consider irregular annulus. 
Using similar conformal transformation, we map the annulus with irregular boundaries from 
(uv) plane onto regular annulus in (xy) plane with inner radius £ and outer radius R = 1. 
Proper conformal transformation looks as follows: 

z + bz 2 + ce t5 z 3 



w(z) 



(A10) 



v /l + 2& 2 (l + £ 2 )+3c 2 (l + £ 2 + £ 4 

For this kind of billiard expansion of ty(r,8) from Eq. (1A5|) is in terms of eigenstates 
bj(r,6) for regular annulus: 

J\l-a\(l\l-a\,nO 



J\l-a\{l\i- a \ >n r) 



N|!-a|(7|J-a|,nO 



N\l-a\(l\l-a\, n r) 



AW 



2ix\l I drr 



J\i- a \{l\i- a \, n r) 



J\l-a\(l\l-a\,nO 
N\l-a\(l\l-a\,nO 



(All) 



iV|i_ a |(7|/_ Q | in r) 
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The counterpart of matrix M for the annulus is: 



(A12) 



+ 5 kh+2 Zce iS l\f I {I + 26 2 (1 + e) + 3c 2 (l + £ 2 + £ 4 )). 



where the integrals l\- 



are defined as: 




drr h+1 



) Y 

lilj 



(A13) 
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